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ABSTRACT 



Context. Studying the composition of dust in the interstellar medium (ISM) is crucial in understanding the cycle of dust in our galaxy. 
Aims. The mid-infrared spectral signature of amorphous silicates, the most abundant dust species in the ISM, is studied in different 
lines-of-sight through the Galactic plane, thus probing different conditions in the ISM. 

Methods. We have analysed 10 spectra from the Spitzer archive, of which 6 lines-of-sight probe diffuse interstellar medium material 
and 4 probe molecular cloud material. The 9.7 /jm silicate absorption features in 7 of these spectra were studied in terms of their shape 
and strength. In addition, the shape of the 18 /jm silicate absorption features in 4 of the diffuse sightline spectra were analysed. 
Results. The 9.7 jim silicate absorption bands in the diffuse sightlines show a strikingly similar band shape. This is also the case for all 
but one of the 18 yum silicate absorption bands observed in diffuse lines-of-sight. The 9.7 yum bands in the 4 molecular sightlines show 
small variations in shape. These modest variations in the band shape are inconsistent with the interpretation of the large variations in 
Tgj/EQ-K) between diffuse and molecular sightlines in terms of silicate grain growth. Instead, we suggest that the large changes in 
Tg -]/E(]-K) must be due to changes in £(J-K). 

Key words. ISM: dust, extinction - ISM: evolution - Techniques: spectroscopic - Infrared: ISM 



1. Introduction 

The composition of dust in the interstellar medium (ISM) is a re- 
sult of a variety of different processes. First of all, dust is formed 
in the circumstellar environments of evolved stars, for example 
. asymptotic giant branch (AGB) stars. Depending on their stage 
' of evolution these stars produce either oxygen-rich (e.g. silicates 
' and oxides) or carbon-rich dust (e.g. amorphous carbon, graphite 
and silicon carbide). The newly-formed Stardust enters the ISM 
through stellar winds or supernova explosions. This material is 
rapidly mixed with other gas and dust in the ISM, where it is pro- 
cessed by a number of mechanisms such as shock waves driven 
by supernova explosions, high energy radiation and high veloc- 
ity collisions amongst grains. Finally, during the formation of 
a new star and planetary system the dust is further processed. 
Therefore, the composition of dust in the ISM is a reflection of 
the constant formation, mixing, processing and destruction of 
different dust species. Studying the composition of interstellar 
dust and its spatial variations throughout the galaxy is crucial 
to get a better understanding of the material from which, for in- 
stance, our solar system was made. 



A good way to study the composition of interstellar dust is by 
means of infrared spectroscopy. The light from the bright back- 
ground star is attenuated by the interstellar dust in front of it. 
This extinction is wavelength dependent and is caused by scat- 
tering as well as absorption by the dust grains. If the intrinsic, 
unreddened spectrum of the background star is known, the dust 
extinction as a function of wavelength can be derived. 

The most abundant dust species in the ISM are amor- 
phous silicates, which cause two prominent abso rption f eatures 
at about 9.7 and 18 A<m (e.g. Stein & Gillett, 1971; Rieke, 
T974tlGillettetal.L [T97l iMcCarthv et al., 1980; Roche, 1989 



Bowev et alJ 



2004; Chiar & Tielens, 2006 



, ■1998HKemperetal , , _ 

Min et al.L 1200 7*). In this article we will discuss both the 9.7 
fim and 18 fim bands. Interstellar silicate dust is mostly a 
mixture of amorphous silicates with an olivine (0/Si=4) or 
pyroxene (0/Si=3) composition (e.g . Davl Il974t Gillett et al.L 



[1971 iKemper et ai.L 125^1 iMin et all l2007[ IChiar et al.L l2007h 



Although olivines and pyroxenes are by definition crystalline sil- 
icates, in this paper we will also use these names for amorphous 
silicates with the same stoichiometric composition. From previ- 
ous studies we know that the composition of interstellar silicates 
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target name 

StRS 136 
StRS 164 

SSTc2d J 182835. 8+0026 1 6 
[SVS76] Ser 9' 

StRS 354 
Elias 3-13^ 

SSTc2dJ163346.2-242753 
G323. 2 103-00.3473'* 

G343.6142-00.1596-* 

G345. 3650-00.40 15* 



other 1 b Program ID spect. diffuse (D) or 9.7 or 18 yum 

designations (degrees) (degrees) type molecular (M) silicate feature 



2MASS J17475608 
-2911439 
2MASS J18161876 
-1635468 

2MASS J 18294508 
+0118469 
IRAS 20273+3740 
IRAS 04303+2609 

2MASS J15285631 
-5653045 
2MASS J16585873 
-4220543 
2MASS J 170706 16 
-4118140 



0.040937 

14.213568 

30.712112 
31.626616 

76.972032 
172.694191 
354.158006 
323.211041 

343.613830 

345.364763 



-0.566892 

-0.002043 

5.280682 
5.423514 



3616 
3616 

139 & 172-179 
139 & 172-179 



0.160012 
-0.401650 



3616 
3616 



B8-A9I- D 
B8-A9I- D 



-0.635686 20294 07~B3^ 

-14.498000 139 & 172-179 K2IlP 

15.606477 139 & 172-179 ... 

-0.347425 3616 



M (Serpens) 
M (Serpens) 

D 

M (Taurus) 



9.7 

9.7 

9.7 
9.7 



9.7 & 18 
9.7 



M (p Ophiuchi) 9.7 

D IB 

D IB 

D IB 



Table 1. Overview of the selected spectra from the Spitzer archive. Listed are the target name and the Galactic longitude and latitude. 
The Spitzer ID of the program that contains the observation and the spectral type of the source is listed and it is indicated if the 
observation was done in a diffuse (D) or molecular (M) sightline. In the case of molecular sightlines, the corr esponding molecula r 
cloud is listed. Finally , we list wheth e r the 9.7 or 18 um siUcate feature was extracted from these spectra. IStrom et alj ( Il976l) . 
^Rawlings et aP (l2QQ0l) . ^ lEliasI (Il978h . lEgan et alJ (l2QQ3h . 



can vary, dep ending on the location in the ISM. For example 
iDemvket ail (1999. 2000, 2001) showed that the dust around 
protostars has a higher O/Si ratio than the newly-formed dust 
observed around AGB stars. 

IChiar et aP (l2007h compared the strength of the 9.7 fim sil- 
icate absorption feature with the near-infrared colour excess 
^(J-K) in both diffuse and molecular sightlinefl The 9.7 fim 
silicate feature, Tgj, is due to the Si-O stretching mode in the sili- 
cates, and its strength depends on the silicate optical depth, while 
the overall extinction in the visible through near-infrared is prob- 
ably due to a combination of carbonaceous material (e.g. amor- 
phous carbon; Draine & Lee, 1984; Co mpiegne et al.^, 2010) and 
iron, the latter either in metallic form dKemper et aU 120021) or 
as a c ationic part of the silicate lattice, the so-called dirty sil- 
icates d Jones & Merrili Il976l) . For diffuse sightlines there is a 
tight Unear co ii'elation between T9.7 and the visual extinction 
(Av) (iRoche & Aitkea 119841: IWhittetL l2003h . A similar corre- 
lation exists between T9.7 and the near-infr ared extinction rep - 
resented by the J-K colour excess, £(J-K) (IChiaret aUlIOOTl) . 
However, for the molecular sightlines the diffu se ISM correla- 
tion fails (IWhittet et ali[T988tlChiar et alil200l . indicating that 
something happens to the dust grains when they enter a molec- 
ular cloud. Molecular sightlines, in general, probe greater inter- 
stellar densities than diffuse sightlines and thus, the observed 
variations are most likely density related. 

An obvious explanation is grain growth via coagulation tak- 
ing place in molecular clouds. Indeed, this causes a decrease 
in the strength of the 9.7 yum silicate feature (see Sect. HI, but 
grain growth also has an effect on the near-infrared extinction 
(i.e. ^(J-K)). Moreover, grain growth alters the shape of th e 
9.7 fim silicate feature drastically (e.g. iBouwman et all l200lh . 
Therefore, an accurate measurement of the band shape of the 9.7 
fim silicate band in diffuse and molecular lines-of-sight can set 



' Diffuse and molecular sightlines are defined as lines-of-sight which 
pass mostly, but not exclusively, through diffuse or molecular cloud ma- 
terial, respectively. 



interesting limits on the importance of grain growth as an expla- 
nation for the large changes in extinction detected so far 

In order to study the spectral signature of interstellar sili- 
cates in both diffuse and molecular sightlines, we have analysed 
10 low -resolution spectra with high S/N ratio, taken with the 
Spitzer space telescope. The high quality of the spectra allows 
an accurate comparison of band shapes along different lines-of- 
sight at a variety of Galactic longitudes. Of these spectra, 7 were 
used to analyse the shape and strength of the 9.7 fim silicate ab- 
sorption feature and 4 were used to analyse the 18 fim silicate 
absorption profile. Only a few earlier studies have analysed the 
interstellar 18 fim band, and always in conjunction with th e 9.7 
fim band (e.g. iDemvk et al.llI999tlChiar & Tielensll2006h . 

This paper is organised as follows: Section 2 describes the 
sample selection and the observations. In Sect. 3 the data anal- 
ysis is described and the results are presented. These results are 
compared with model calculations in Sect. 4, and in Sect. 5 we 
discuss the implications of our results for grain processing in 
molecular clouds. Finally, the conclusions are summarised in 
Sect. 6. 

2. Observations 

2.1. Sample selection 

Our sample consists of a selection of low -resolution IRS 
(InfraRed Spectrograph; i Houck et aI.L |2004|) Spitzer spectra 
from the Spitzer programs with Program IDs: 139, 172-179, 
3616, 20294, 3336, 20604. These programs have, amongst other 
things, observed background sources with large amounts of fore- 
ground dust and were thus relevant for this study. The back- 
ground nature of the dense cloud sources was based on their 
placement in colour-colour diagrams which separated out stars 
reddened by interstellar, rather than circumstellar dust. For the 
diffuse ISM lines-of-sight, the background nature is confirmed 
via their strong reddening an d optical spectra that c onfirm their 
spectral and luminosity class (iRawlings et al.Ll2000l) . 
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A first selection was made based on the presence and 
strength of the 9.7 yum and 18 /im siHcate absorption features and 
the signal to noise ratio (SNR) of the spectra. From the result- 
ing list 10 spectra were selected keeping in mind the following 
science goals: (1) Investigating the optical depth at 9.7 fim, rg 7, 
versus the infrared colour excess, ^(J-K), (2) Comparing the 
shape of the 9.7 jum silicate features through different sightlines, 
(3) Comparing the shape of the 18 jim silicate features through 
different sightlines. 

To reach the first two science goals the 9.7 fim silicate fea- 
ture needs to be extracted from the spectrum. However, many of 
the spectra of background stars in the Spitzer archive show an 
SiO gas phase fundamental stretch band from about 8 to 1 1 yum 
coming from the intrinsic spectrum of the observed background 
star (see Fig.[T]). This photospheric gas phase SiO band is typical 
for late type K and M giants and because these are intrinsically 
bright stars, they are good candidates for background sources. 




photosphere 

80 (SiO) < < < < 

5 10 15 20 25 30 

Fig. 1. Spectrum of an M-giant taken with the Spitzer Space 
Telescope (Chiar, unpublished Spitzer data). The extinction to- 
wards this star is Ay ~ 10 mag. Plotted is the observed flux (in 
Jy) of the star multiplied by the wavelength (in fim) squared 
(Fy X A^) versus the wavelength (A). This choice implies that a 
horizontal line is expected for an unreddened photosphere. The 
two most prominent absorption bands at 9.7 and 18 fim are due 
to interstellar silicate absorption. The smaller band at 8 jum is 
caused by gas phase SiO in the photosphere of the M-giant itself 
and clearly affects the 9.7 fim silicate feature, thus hampering 
the detailed analysis of the shape of the interstellar silicate band. 

Since the shape of the SiO band depends on the spectral 
type of the background source, which in most cases is not 
known, accurately coiTecting for this band is not easy (see also 
Sect. 13.3.2b . The profile of the SiO fundamental mode is highly 
sensitive to the temperature and the structure of the stellar at- 
mosphere, re quiring modelling beyond that of a classical atmo- 
sphere (e.g. iTsuii et all Il994t) . Therefore, 7 spectra were se- 
lected that show no or very little SiO absorption. Of these 7 
spectra 3 are diffuse sightlines and 4 are molecular sightlines. 
All but one of these spectra were not suited to analyse the 18 fim 
silicate absorption feature, either because the spectrum did not 
span the full spectral range or the quality of the long wavelength 
part of the spectrum was not high enough to reliably extract the 
18 fim feature. Therefore, we have selected 3 more spectra from 
the Spitzer archive to analyse the 18 /im silicate absorption fea- 
ture, which in turn were not suited for the analysis of the 9.7 fim 
feature, because they show strong photospheric gas phase SiO 



bands (see Fig. [U. Unfortunately, we did not find any molec- 
ular sightline spectra in the Spitzer archive that span the long 
wavelength part of the Spitzer spectral range and therefore, we 
could only analyse the 18 fim silicate absorption feature in 4 dif- 
fuse sightlines. An overview of the selected spectra is given in 
Tabled 

2.2. Data reduction 

Basic calibrated data (BCD) were obtained from the selected ob- 
servations using the Spitzer data reduction pipeline version S15 
and post-BCD ID spectra were created by carrying out back- 
ground subtraction. The resulting spectra are shown in Fig. |2] 
and|3] We plot these as (Fy x A^) versus A so that an unreddened 
photosphere without photospheric absorption bands shows up as 
a horizontal line, making it easier to recognise the effect of ex- 
tinction on the spectral slope. 

2.3. The spectra 

The 9.7 fim silicate feature is present in all spectra (Figs. |2] and 
|4]l. The photospheric gas phase SiO band is not apparent in the 
spectra, but for only three sources we can be sure that it is indeed 
not present, since they are early type supergiants (see Table [TJ. 
For the other sources, however, a very weak SiO band might still 
be present. 

Ice features are reported to b e absent from the 
Spitzer spectra of diffuse sightlines (Ivan DishoeckL 120041 : 
iBoogert & Ehrenfreund, 2004). Water ice, for instance, has a 
feature at about 6 fim and CH^O H ice together with N H| ice 
has a feature at about 6.85 fim (e.g. I Boogert et al.Ll2008l) . which 
we do not detect in the diffuse sightline spectra in this sample. 
This is as expected, since the conditions in the diffuse ISM 
are not suitable for the formation of ices. Ice features are also 
not detected in the spectra towards sightlines with Ay < 10, 
although the 6.0 fim water ice resonance appears only as a weak 
absorption towards sightlines with higher column densities 
(Ay > 10). Therefore, the 9.7 fim silicate features are probably 
also not strongly affected by ice bands, such as the 12 fim water 
ice librational band. Our sightlines do probably not sample the 
densest parts of the molecular clouds, where ices are expected 
to be abundant. 

The spectra in Fig.|3]show the 18 fim silicate feature. In this 
spectral range the signal to noise level of the data is rapidly de- 
creasing due to the lower flux levels. We estimate the S/N of the 
spectra to be in the range 15 to 25. This will affect the accuracy 
of the extracted band profiles. 

3. Data analysis and results 

3.1. Extraction of the 9. 7 and 18 fim silicate absorption 
features 

The total extinction by interstellar dust can be divided into 
a "continuum" component (i.e. smoothly varying with wave- 
length), and the silicate absorption features, the shape and 
strength of which we wish to extract. To do this a continuum 
needs to be specified, which represents the intrinsic spectrum 
of the background star reddened by foreground dust. We as- 
sume the intrinsic spectrum of the background star to be approx- 
imately a Rayleigh- Jeans tail, described by Fy oc A^^, while the 
featureless dust extinction can be approximated by a power-law 
with wavelength as well. The resulting continuum, against which 
the spectral features are seen in absorption, is thus represented 
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10 
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Fig. 2. Flux multiplied by wavelength squared spectra of the sightlines used for the 9.7 fim silicate feature study together with the 
continuum (see main text). 
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X (iim) 



160 
140 
120 - 
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G343.6142-00.1506 
(diffuse sigfit line) 



14 



20 22 24 26 

Xitim) 



120 
115 
110 
105 
100 
95 
90 
85 



150 
140 
130 
120 
110 
100 



G323.21 03-00.3473 
(diffuse sigfit line) 



20 22 

Xitim) 



90 - 
80 



G345.3650-00.4015 
(diffuse sigfit line) 



20 22 24 26 

Xitim) 



Fig. 3. Flux multiplied by wavelength squared spectra of the sightlines used for the 18 jim silicate feature study together with the 
continuum (see main text). Error bars: 2, 3, 3 and 5 JyXyum^ for StRS354, G323, G343, and G345 respectively. 



by a power-law. This approximation is good enough for our pur- 
poses, since the widths of 9.7 jum and 18 fj.m silicate absorption 
features each cover a relatively small wavelength range (~4 and 
~13 jum, respectively). 

The continuum for the 9.7 fim silicate feature is defined by 
the flux levels at 7.5 and 12.3 fj.m, where the first continuum 
point (at 7.5 jum) is chosen because it does not overlap with po- 
tential photospheric gas phase SiO bands that might be present 
in the spectra. The resulting slopes are listed in the last column 
of Table|2] The continuum for the 1 8 fim silicate feature is deter- 



mined separately and defined by the values at 14.2 and 27.0 jum, 
thus making sure both continuum points are within the wave- 
length range of the long wavelength module of Spitzer. 

Subsequently, the optical depth profiles are calculated using: 

r = -log(|^) (1) 

^ ^ v,cont ' 

where Fv obs is the observed flux and /^v,cont is the flux of the con- 
tinuum. The resulting optical depth profiles are shown in Figs.|4] 
and|5] and the peak values for the 9.7 yum feature are also listed 
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9.5 10 10.5 11 11.5 12 

X (Hm) 




;.5 9 9.5 10 10.5 11 11.5 12 



X (|J,m) 




5 9 9.5 10 10.5 11 11.5 12 



X (|im) 



Fig. 4. The extracted 9.7 fim optical depth profiles. 



in Table[2] A weak correlation exists between T9.7 and the power- 
law index of the fitted continuum, with a coiTelation coefficient 
of -0.59. 

The choice of the continuum is important because it af- 
fects the shape of the extracted band profile and the subsequent 
analysis. We chose to use a systematic approach by always us- 
ing the same wavelengths at which the continuum is defined. 
We note that the three sources with the highest extinction in 
our sample (all are molecular sightlines) show such a sharp 
drop in flux shortward of 8 jim that this falls below the ex- 
trapolated continuum. We have investigated the impact of this 
choice on our analysis by choosing a different continuum for 
two objects, StRs 164 and SSTc2dJ163346.2-242753 in the 
9.7 fim region. The strength of the band changes by at most 
10%, when normalised to unity at a wavelength of 10.5 fim, 
while the shape remains virtually unaffected. This is shown for 
SSTc2dJ163346.2-242753 in Fig.|6] We conclude that this dif- 
ferent choice of continuum has a minor effect on the outcome of 
our analysis. 



3.2. T9.7 versus E(J-K) relationship 

IChiar et aH (l2007h investigated the relationship between the 
strength of the 9.7 fj.m silicate feature, T97, and the near- 
infrared colour excess, ^(J-K). For diffuse sightlines there is 
a tight linear coiTelation between these two parameters with 
T9.7/£(J-K)=0.34 (Roche & Aitken, 1984; Whittet, 2 00l, but 
for the molecular s i ghtUn es this correlation fails dChiar et al.L 
I2007t IWhittet etallll988h . The extinction law within the NIR 



wavelength regime appe ars to remain unchanged w ith increas- 
ing column density (e.g. iRoman-Zufliga et all l2007h . 

To compare the sources in our sample to these previous re- 
sults, T9.7 was determined from the 9.7 fim silicate features in 
our sample. An estimation of the error on these measurements 
was done by taking different reasonable continua and determin- 
ing the effect on T9 7. 

The ^(J-K) colour excesses were determined b y taking the 
(J-K) colour observed by 2MASS (Skrutskie et al., 2006) and 
subtracting the intrinsic (J-K) colour from Koornneefi(J983l) . If 
the spectral type of the source is unknown, we assume it is a GO- 
M4 giant, with an average intrinsic (J-K) colour of about 0.81 
magnitudes (consistent with Chiar et al. (2007)). The errors were 
estimated by taking into account the uncertainty in the spectral 
type and thus the intrinsic (J-K) colour. The results are listed in 
Table|2]and Fig.Qshows these results toge ther with the Gal actic 
Centre line-of-sight and the results from Chiar et aTl (120071) for 
the diffuse and molecular sightlines separately. 

Our measuremen ts are in agreement with the results found 
bv lChiar et al ] (I2007h as can be seen in Fig. |7] Please note, that 
two of our measurements are for for sightlines also in the sample 
oflChiar et al. (2007), namely: EHas 3-13 and [SVS76] Ser9 and 
we have replaced their measurements by ours in Fig.|2] The two 
measurements for [SVS76] Ser 9 agree within th e uncertainty. 
For Elias 3-13 the measurement was done by IWhittet et al.l 
(il988D from a relatively noisy, ground based spectrum, so the 
error on their measured T9 7 is quite large (approximately 0.08 
in r). Moreover, they did not correct for photospheric gas phase 
SiO. 
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target name 


Spectral type 


Diffuse/Molecular 


(J-K),,,,' 


(J-K)oftj 


£(J-K) 




power-law index 


StRS 136 


B8-A9I2 


D 


0.12 


1.78 


1.66±0.1 


0.65+0.05 


-2.00 


StRS 164 


B8-A9I- 


D 


0.12 


1.78 


1.66±0.1 


0.62+0.05 


-2.03 


SSTc2dJ182835. 8+002616 




M (Serpens) 


0.81 


3.78 


2.97±0.40 


0.73+0.05 


-2.25 


[SVS76] Ser 9 




M (Serpens) 


0.81 


4.33 


3.52+0.40 


0.74+0.05 


-2.27 


StRS 354 


07-B32 


D 


-0.09 


1.88 


1.97+0.1 


0.58+0.04 


-1.73 


Ellas 3-13 


K2III3 


M (Taurus) 


0.73 


3.01 


2.28+0.05 


0.42+0.04 


-2.15 


SSTc2d J 1 63346.2-242753 




M (p Ophiuchi) 


0.81 


2.45 


1.64+0.25 


0.45+0.05 


-1.28 



Table 2. The results for rg 7 and ^(J-K) measured for the sources in our sample. The target name, spectral type (if known) and 
the type of sightline (diffuse, D, or molecular, M) are listed. If the spectral type is unknown we assume it is a G0-M4 giant. The 
corresponding intrinsic J-K colour, (J-K),„,, and the J-K colour as observed by 2MASS, {i-K)obs, are also listed. The difference 
between these two values gives the colour excess ^(J-K). The last two columns list the the r neasur e d optical depth at 9.7 ju m , T9.7, 
and th e index of the power-law used to determine the continuum around 10 jjm. iKoornne^ (Il983l) . ^Rawlings et al.l (|2000|) . feliasl 
(fT978h ) 



We also included the Galactic Centre l ine-of-sight iri Fig. [71 
(upper panel). F or this we used rg 7=3.25 ('Kemper et al.', '20()4t 
iMin et aUl2007l) and T9.7/£'(J-K)=0.7 (Hoche & Aitken, 1985^. 
The Galactic Centre line-of-sight is exceptional, since it is the 
only one for which T97 versus ^(J-K) lies above the diffuse 
ISM correlation. This sightline probes the local ISM as well as 
the inner parts of our Galaxy in contrast to the other sightlines 
which only probe the local ISM. Since there are relatively fewer 
C-rich stars (as compared to O-rich stars) around the Galactic 
Centre than in the outer regions of our Galaxy, the relative abun- 
dance of carbo n-rich dust is probably also lower around the 
Galactic Centre dRoche & Aitkenl[l98 5'). This would explain the 
relatively high T9 7/£'(J-K) ratio for this particular line-of-sight, 
since carbon rich dust is a major component of the near-infrared 
extinction. 

Furthermore, there is one molecular cloud source (indicated 
by the arrow in Fig. |7]) that seems to agree very well with the 
diffuse ISM correlation despite the high value for ^(J-K), while 
other molecular sightlines with ^(J-K) > 2 deviate from this 
correlation. This is SSTc2dJ182852.7-H002824, w hich was ob- 
served with Spitzer in the c2d legacy program (IE vans et al.L 



|2003|) . Its reduced spectrum is shown in Fig. [8] Prominent ice 
features are present at about 6 and 6 .85 um caused by H7O 
and CH3OH -I- NHt ice r espectively (iBoogert & Ehrenfreundl 



I2004tlvan Dishoeckil2004 . The latter also has a feature at about 
9.7 /im, which increases the depth of the 9.7 jiva silicate fea- 
ture leadin g to an overestima te of the absorption due to silicates. 
Therefore. IChiar et alj (12007 ) apparently overestimated the ab- 
sorption due to siUcates. Correcting for ice absorption would 
move the T9 7 point downward in Fig. |7] which is in agreement 
with the results for the other dense cloud sources. We note the 
large differences in the strength of the ice bands between the 
spectrum in Fig. [8] and the 4 molecular cloud sightlines shown 
in Fig.|2] for similar values of ^(J-K). Apparently, the observed 
changes in the 79 7 versus ^(J-K) relationship are not related to 
ice formation. 



3.3. Comparison of the 9.7 urn silicate profiles 

The shape of the 9.7 /urn silicate absorption features in our sam- 
ple (see Fig. |4]i will be compared and discussed for the diffuse 
and molecular sightlines separately below. Throughout this arti- 
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Fig. 6. The effect of choosing a different continuum on the shape 
of the 9.7 fim siHcate band for the molecular sightline source 
SSTc2dJ1633462-242753. The rising continuum in the 18 yum 
wavelength range is due to the "continuum extinction" which 
decreases with increasing wavelength, upper, two choices of the 
continuum normalisation; lower: comparison of the old and new 
band shape. 



cle we will use the fe ature determined for the Galactic Centre by 
iKemper et al.l (l2004l) as a reference. 



3.3.1. Diffuse sightlines 

Figure |9] shows the extracted features for the diffuse sightlines 
in our sample toget her with the Galactic Centre feature from 
iKemper et al.l (12004 '). All features are normalised to unity at 10.5 
fim, because the contribution of any photospheric gas phase SiO 
bands that might be present in the spectra is negligible from 
about 10.5 i-im onwards. This effect is not important for the dif- 
fuse sightlines, since the background sources are known to be 
early type supergiants without photospheric gas phase SiO, but 
it will matter for the molecular sightlines which we will discuss 
later. 

From Fig. |9]we conclude that the profile of the 9.7 yum sil- 
icate features in different diffuse sightlines is very similar. In 
Fig. [To] the difference between our observed silicate features 
and the Galactic Centre feature (Kemper et al., 2004) is shown. 
For both StRS 164 and StRS 354 there are no strong varia- 
tions present that exceed the noise. Only StRS 136 shows a 
small difference with the Galactic Centre feature between about 
8 and 9.7 /im, but this never exceeds 0.15 in the normalised 
T. The two small absorption features at about 8.7 and 11.3 
yum in Fig. [10] (upper panel) cannot be trusted, because PAHs 
(Polycyclic Aromatic Hydrocarbons) have prominent features at 



Diffuse ISM 
sight lines 




Chiar et al. 2007, and references therein 
new data 
diffuse iSM correiation 



Molecular cloud 
sight lines 




Fig. 7. The measured optical depth at 9.7 yum, T9 7, as a func- 
tion of the J-K colour excess, ^(J-K), for diffuse sightlines 
(upper panel) and molecular sightlines (lower panel). The dots 
represent data from Chiar et al. (2007, and references therein) 
and the crosses represent our measurements. The dashed line is 
the coiTelation between Tgj and ^(J-K) found for the diffuse 
ISM (Chiar etal., 2007). The source indicated by the arrow is 
SSTc2dJ182852.7-i-002824. This source is remarkable, since it 
is the only source with £'(J-K)>2 that is still in agreement with 
the diffuse ISM correlation, probably due to the effects of ices 
on the estimated band strength. 



those wavelengths and the background subtraction probably was 
not accurate enough in this case. 

In conclusion, there are only small differences in the shape 
of the 9.7 pm silicate feature in the diffuse sightlines and the 
Galactic Centre feature, even though these sightlines are sepa- 
rated by up to 77 degrees in Galactic longitude (see Table]!). 



3.3.2. Molecular siglitlines 

Figure ]TT] shows the extracted features for the molecular sight- 
lines in our sample. All features are normalised to unity at 10.5 
pm as was done for the diffuse sightlines. From an eyeball com- 
parison of the 9.7 pm silicate features in Fig.fTTIwe conclude that 
these features in molecular sightlines differ substantially from 



8 



J.M. van Breemen et al.: Silicate absorption profiles in molecular and diffuse interstellar sightlines 





SSTc2d_J 182852.7+002824 —I — 


/ CH30H \ 




-' + NH4*? 




H20 






NH3? 




CH30H 



7 
6 

X 3 
2 
1 



5 6 7 8 9 10 11 12 13 14 15 
'k (urn) 

Fig. 8. The spectrum observed towards 
SSTc2dJ 182852.7+002824. Absorption features caused 
by ices on top of the silicate feature are indicated. The absorp- 
tion feature caused by CH3OH ice at about 9.7 //m makes the 
silicate feature seem deeper than it really is. Consequently, the 
strength of the 9.7 jum silicate feature is overestimated. 
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Fig. 9. The 9.7 jim silicate features in terms of optical depth for 
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malised to unity at 10.5 jjm. 
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the absorption feature observed towards the Galactic Centre. The 
slopes of the long wavelength side of the features (i.e. A > 10.5 
fj.m) are the same as that of the Galactic Centre. If we assume 
that the profiles are indeed the same from 10.5 yum onward, then 
the molecular sightlines show excess absorption between 7.5 and 
10.5 yum. Alternatively, if we were to normalise the 9.7 yum sil- 
icate features at their peak, the difference between them can be 
interpreted as a shift towards shorter wavelengths as compared 
to the diffuse 9.7 jim silicate features. Also, the long wavelength 
part of the features show less overlap than when they are nor- 
malised at 10.5 yum. 

In Fig. [12] the difference is plotted between the observed 
features in the molecular sightlines in our sample and the 
Galactic Centre feature. The observed difference between the 
9.7 yum silicate profiles in molecular sightlines and the Galactic 
Centre sightline can be caused by three effects: (i) The presence 
of a photospheric gas phase SiO band in the spectrum of the 
background star, (ii) The presence of ices along the line-of-sight. 
(iii) A change in the dust properties of the amorphous silicates. 
We will discuss these three effects separately below. 

(i) Photospheric gas phase SiO band: Since the spectral 



Fig. 11. The 9.7 yum silicate features in terms of optical depth 
for the molecular sightlines together with the fe ature observed 
towards the Galactic Centre dKemper et all I2004D . The features 
are normalised to unity at 10.5 fim. 



type for 3 of the 4 sources is unknown, there might be a small 
photospheric gas phase SiO band in the spectra. SiO is the 
most probable source of photospheric spectral structure in 
cool stars in the 8 yum wavelength range. Indeed, Elias 3-13 is 
known to be a K2 giant and, hence, the spectrum of this source 
will have a weak SiO band. Careful inspection of the 9.7 yum 
features in Fig. |4] shows that the photospheric gas phase SiO 
band is indeed present in all but one of the molecular sightlines 
as a blue shoulder between 7.5 and about 8.5 yum. Only for 
SSTc2dJ182835.8-i-002616 the SiO band seems to be absent, 
since r is zero at about 8 yum, where the SiO band peaks. So for 
this source, the observed difference with the Galactic Centre 
feature cannot be explained by the presence of a gas phase SiO 
band in the intrinsic spectrum of the background star. 

However, Fig. [12] shows that the difference in the shape of 
the 9.7 yum silicate features, in most cases, cannot be explained 
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Fig. 12. The difference between the observed features in the 
molecular sightlines in our sample and the Galactic Centre fea- 
ture derived by Kemper et al, (2004 ). (The Elias 3-13 spectrum 
has not been corrected for the photospheric gas phase SiO band 
in this figure) 



only by SiO. This is because gas phase SiO normally peaks at 
about 8 yum, whereas Fig. [12] shows features that peak at about 
9-9.5 yum, except for SSTc2dJ163346.2-242753. 

Since the spectral type of Elias 3-13 is known we were 
able to extract the photospheric gas phase SiO band for this 
source using a mod el atmosphere spectrum calculated by one 
of us (iDecin et al.L [19 97. 2004) based on MARCS models 
jGustafsson et al.L 1 19751 1200 8) (see Fig.[T3] upper panel). 

The photospheric gas phase SiO band in terms of optical 
depth was extracted from this model spectrum by taking a 
power-law continuum from 7.4 to 12.3 fim and using equation 
1 (see Fig. [13] lower panel). The resulting optical depth profile 
is mostly due to SiO, but a contribution from OH is present 
at /I > 9 /im, and OH absorption dominates at > 10 /urn. 
Subsequently, we subtracted this optical depth profile from the 
observed 9.7 fim feature. The correction never exceeds 0.03 
in T. The result is shown in Fig. [14] Correcting for SiO in 
this spectrum is not sufficient to explain the difference in the 
shape of the 9.7 yum silicate feature between Elias 3-13 and the 
Galactic Centre. 



(ii) Ices: Different types of ices could affect the shape 
of the 9.7 ^m sihca t e fea t ure in mole c ular sightlines 
dBooge rt & Ehrenfreundi l2004t Ivan DishoecH |2004 . Indeed, 
[McClurJloOsF proposed ice growth as cause for the change 
in extinction law between diffuse and molecular sightlines. For 
instance, CH3OH ice has a strong feature at about 9.7 jjm (see 
also Fig. [8]). Since the spectra show no strong absorption bands 
at about 6 //m and 6.85 yum, caused by water ice and CH3OH 
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Fig. 14. The 9.7 //m silicate feature observed towards Elias 3- 
13 (dashed line) together with this feature corrected for pho- 
tospheric gas phase Si O (dotted line) using a m odel spectrum 
calculated by L Decin (iDecin et all [19971 12004 . The Galactic 
Centre feature (iKemper et al.L 12004 1)) is also plotted (solid line) 
and is representative for diffuse sightlines. 



ice together with NH4 ice respectively (see Fig. [2]), we do not 
expect the 9.7 fim. silicate features to be affected much by ice 
features. 

For the same reason we do not expect a strong H2O libration 
band at about 12 yum (Boogert et al., 2008). Moreover, the long 
wavelength side of the 9.7 jjm silicate features (from about 10.5 
jim onward) in our molecular sightlines is very similar to that 
of the diffuse sightlines indicating that there is indeed no H2O 
libration band present in any of the spectra. 

In Fig.[l2] only SSTc2dJ182835.8-H002616 shows a feature 
that peaks at about 9.7 /^m, that may be compatible with the 
CH3OH ice band. Figure [15] shows the difference between the 
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SSTc2dJ182835.8+002616 and the Galactic Centre feature de- 
rived by .Kemper et al.. (.2004.) (solid line) together with the ab- 
sorption profile of a mixture of CH3OH and CH3 ice (dashed 
line). 



9.7 jum feature of SSTc2dJ182835. 8+002616 and the Galactic 
Centre feature together with the absorption profile of a mixture 
of CH3OH and NH3 ice. In this case the observed difference 
could be explained by the presence of ices. On the other hand 
the spectrum of SSTc2dJ182835.8+002616 does not show an 
absorption band at 6.85 /im, caused by CH3OH ice together with 
NH+ ice. 

For Elias 3-13, however, the difference between its 9.7 fim 
feature and the Galactic Centre feature peaks at about 9 yum even 
if we correct for the photospheric gas phase SiO band. Moreover, 
it is a relatively broad feature covering about 2-3 fim in wave- 
length. There are, to the best of our knowledge, no ice species 
that produce such an absorption profile. 

(iii) Dust properties of the amorphous silicates: The shape 
of the 9.7 jum silicate absorption feature depends on the dust 
properties of the amorphous silicates, such as the size, shape 
and chemical composition of the dust grains. The observed in- 
crease in absorption between 7.5 and 10.5 fim (i.e the short 
wavelength side of the feature) could, for instance, be the result 
of a decrease in the O/Si ratio in the amorph ous silicates, go- 
ing from an olivine to pyroxene composition (Ossenkop f et al.L 
Ll992; Speck et al., 2008). In Sect.|4]we will return to this point. 

From the discussion above we conclude that it is not likely 
that gas-phase SiO absorption or ice absorption are the dominant 
cause of the difference in the band shape of the 9.7 fim silicate 
absorption band in diffuse and molecular lines-of-sight. Changes 
in the dust composition may be important in explaining the dif- 
ference. 

3.4. Comparison of the 18 fim silicate profiles 

The 18 fim silicate absorption profiles are extracted for 4 diffuse 
sightlines (see Sect. 13. lb . Because we could analyse the 9.7 fim 
feature in only one of these sightlines, it is not possible for us to 
use the feature strength ratio between the two features as a tool 
to put con straints on grain prope rties, such as composition and 
size (e.g. 'O ssenkopf et al.LTl992 ). Figure [16] shows the 18 ;um 
profiles together with a modelled Galactic Centre 18 fim profile, 
based on the grai n properties de rived from the 9.7 fim fit in the 
same sightline ( Min et al.Ll2007h . The features are normalised to 
unity at their peak. 

The 18 fim silicate profile of G345.3650-00.4015 is the only 
one that does not agree with the extrapolated Galactic Centre fit 
within the noise and we have found no reasonable explanation 
for this particular profile being so different. For this sightline, we 
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Fig. 16. The 18 fim silicate absorption features observed in 4 
diffuse sightlines together with the modelled profile based o n 
the 9.7 fim fit of the Galactic Centre sighdine (iMin et al.Ll2007l) . 
The features are normalised at their peak. Noise levels vary with 
wavelength and are typically 0.03 in t at the shortest wave- 
lengths, and increase to 0. 10 in r for the longest wavelength. The 
anomalous feature seen towards G345. 3650-00. 4015 is plotted 
with a thicker line than the features in the other sightlines. 



considered a range of olivine and pyroxene silicates with differ- 
ent Mg/Fe ratios, as well as Na- and Al-containing silicates. For 
all these species we considered spherical and irregularly-shaped 
particles, none of which reproduce the required blueness of this 
anomalous feature. Although we suspect that a compositional 
effect causes the shift in wavelength for this feature, the carrier 
remains unidentified. However, the 18 fim features in the other 
sightlines are very similar to th at of the extrapolated Galactic 
Centre fit and iMin et all (l2007h show that this profile is also 
consistent with the 18 fim silicate absorption feature observed 
towards WR 98a. So, except for one source, the 18 fim silicate 
absorption features in different diffuse sightlines are very similar 
and the extrapolated Galactic Centre fit seems representative for 
the diffuse ISM. 

3.5. Summary of observed trends 

For lines-of-sight passing through the diffuse ISM, there is a 
tight linear coiTelation between the strength of the 9.7 fim sil- 
icate absorption feature, rg 7, and the near-infrared colour ex- 
cess, ^(J-K). Moreover, the shape of the observed 9.7 fim sil- 
icate features in different diffuse sightlines and the Galactic 
Centre feature are strikingly similar (see Fig. |9]l. Only StRS 
136 shows a small difference with the Galactic Centre feature 
between about 8 and 9.7 fim of up to 0.13 in r. It could be 
that this line-of-sight not only probes diffuse material but also 
a small amount of molecular material. If this is indeed the case, 
then the Galactic Centre feature is representative for the diffuse 
ISM even th ough this feature wa s extracted using a very differ- 
ent method d Kemper et all l2004l) than applied here. This is re- 
markable, since the Galactic Centre line-of-sight not only passes 
through diffuse material, but also some molecular cloud material 
( Roche. 1988; Whittet, ,2003|). 

Our analysis of the 18 fim silicate absorption feature shows 
the same similarity between diffuse sightlines and the Galactic 
Centre. Except for one source (i.e. MSXLOS17) the 18 fim pro- 
files are consiste nt with the extrapolated Galactic Centre fit from 
lMinetalJ(l2007h . 
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Muss fmction (%) 


Amorphous olivine 


MgFeSi04 


13.8 




Mg2Si04 


38.3 


Amorphous pyroxene 


MgSiO^ 


42.9 


Amorphous Na/Al pyroxene 


NaAlSi206 


1.8 


Crystalline olivine (Forsterite) 


Mg2Si04 


0.6 


Other 


Crystalline silicon carbide 


SiC 


2.6 



Table 3. The chemical composition in terms of mass fraction 
in percent for the best fit model of the Galactic Centre feature 
bv lMin et al.l (l2007h . This model uses a distribution of hollow 
spheres (DHS) with f^ax - 0.7 and an MRN size distribution 
with index -3.5 between a grain size of 0.005 and 0.25 /im. 



We conclude that the diffuse ISM silicates probed by the 
diffuse lines-of-sight studied in this paper and those sampled 
towards the Galactic Centre are remarkably similar. Only the 
T9.7/£'(J-K) of the Galactic Centre sightline deviates from the 
general behaviour. 

A different picture appears when looking at lines-of-sight 
passing through molecular cloud material. The rg? versus 
E(}-YL) relationship breaks down drastically and we observe 
small variations in the shape of the 9.7 jiva silicate absorption 
features. These small changes may be attributed to differences 
in dust properties, but other effects may also play a role (see 
above). Nevertheless, the observations show that the dust prop- 
erties in molecular clouds must be different from those in the 
diffuse ISM. 

To determine whether the change in the rg 7 versus ^(J-K) 
relationship is correlated to the change in the shape of the 9.7 
//m silicate feature, we have measured the 9.7 ^va silicate band 
shape of three sources in Fig.|2]with a wide range in rg 7/£'(J-K). 
However, we found that better quality data than available to us 
in this study are needed to come to firm conclusions about such 
a possible relation. 

In the following we will consider possible explanations for 
the observed variations in both the T9.7 versus ^(J-K) diagram 
and the shape of the 9.7 /im silicate profile. We model the T9 7 
versus ^(J-K) relationship assuming a range of dust properties; 
then we confront these models with the observations. 



4. Model calculations 

To investigate how the T9.7 versus ^(J-K) relationship and the 
9.7 //m silicate absorption profile depend on different dust prop- 
erties we start with the best fit found by Min et al. (2007) for 
the G alactic Centre feature. This fit uses an MRN (Mathis et al., , 
Il977h size distribution (i.e. n{a) oc a ), between grain sizes of 
0.005 and 0.25 jum. The shape of the particles is represented by a 
distribution of hollow spheres (DHS), where /^ax represents the 
size of the cavities in the grains and is a measure for the irregular- 
ity of the particles, and a large r fmax indic ates that the grains are 
more porous and/or irreg ular jMin et all |2005). A DHS model 
with /max = 1 simulates highly irregular particles, while a DHS 
model with /max =0.1 simulates alm ost spherical particles. The 
best fit model from Mi n et al. I (l200l uses /max - 0.7. The chem- 
ical composition of this best fit model is given in Table[3] 



By changing the parameters in this model we are able to in- 
vestigate the effects of grain size, grain shape and chemical com- 
position of different types of silicates and silicon carbide on the 
9.7 fxra silicate absorption profiles. 

To model the relationship between T9 7 and ^(J-K) we have 
included a fraction of graphite (C) in such a wa y that for the 
Galactic Centre the diffuse ISM coiTelation dChiar et all |2007|) 
is obtained. Note that this correlation is not valid for the Galactic 
Centre sightline in particular (see Fig.|7]i, but the 9.7 yum silicate 
absorption profile is representative for the diffuse ISM and thus 
we will consider it as a typical diffuse sightline. A mass fraction 
of about 3 1 % of graphite was needed to meet this requirement. 
Assuming a solar abundance of carbon, this requires about 65 
% of the carbon atoms to be in the form of graphite, which is a 
relatively high percentage. This number should be treated as an 
upper limit; when metallic iron or dirty silicates are present in 
the diffuse ISM, they will contribute to the NIR opacity, and thus 
decrease the contribution from graphite. 

To verify if this model also reproduces the standard ex- 
tinction law at optical wavelengths the coiTesponding value 
for Ry=Av/£(B-V) was calculated. The result was 3.28, 
which is close to the value for the diffuse I SM of about 3.1 
(Sc hultz & WiemeillI975l:ISneden et aUII978l) . We consider the 
following dust properties and investigate their effect on the T9 7 
versus ^(J-K) relationship: (1) The size distribution of the dust 
grains, (2) The shape distribution of the dust grains, (3) The 
chemical composition of the dust grains. In the following we 
will consider graphites and silicates as separate grains and we 
will assume that they have the same size and shape distribution. 

4.1. Grain size 

Due to the high density, dust grains are likely to grow in molec- 
ular clouds. In order to model this, we increase the maxi- 
mum grain size rmax while keeping rmin constant, and calculate 
T9 7/£'(J-K) as well as the shape of the 9.7 /vm silicate feature. 
The result is shown in Fig.[T7l At first, when the maximum grain 
size is increased, T9 7/£'(J-K) decreases down to about 0.28, but, 
as Tmax excccds ~0.6 //m, it starts increasing. This behaviour 
is dominated by the change in ^(J-K) due to growth of the 
graphite grains. With this model, however, we are not able to ex- 
plain the observed flattening of the rg 7 versus ^(J-K) relation- 
ship, because our model curves do not reach the lower right part 
of the Tg 7, ^(J-K) diagram, where the most reddened molecular 
lines-of-sight are located (see Fig.|7]). 

The shape of the 9.7 //m silicate feature does not change 
much as the maximum grain size is increased up to I.O //m. 
When rmax exceeds this value the 9.7 yum silicate profile weak- 
ens and shifts towards longer wavelengths. This is in contrast 
with our observations which show slight excess absorption on 
the short wavelength side of the feature (see Fig.fTTTi. So, if grain 
growth is indeed taking place, it is limited to a grain size of about 
1.0 ;um, and cannot explain the observed variations. 

4.2. Grain siiape 

Figure [18] (upper panel) shows how the Tg 7 versus ^(J-K) rela- 
tionship changes as a function of fmax- With increased sphericity, 
the slope of the relation decreases towards the value observed in 
molecular sightlines (see Fig. |7ll. The effect is small, however, 
and certainly not strong enough to explain the observed break 
down of this relationship in the observations. 
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Fig. 17. upper: The effect of increasing the maximum grain size 
on the the T9 7 versus ^(J-K) relationship. The minimum and 
maximum grains size are given in fims. lower: The effect of in- 
creasing the maximum grain size on the shape of the 9.7 jjm 
silicate feature. The left panel shows the calculated 9.7 fim sil- 
icate absorption profiles in terms of the mass extinction coeffi- 
cient and the right panel shows the same profiles normalised to 
1 at 10.5 fim as was done with the observations. The solid line in 
both panels rep r esent the best fit for the Galactic Centre feature 
from (Minetal](l2007h . 



The lower panel of Fig. [18] shows the effect of changing the 
grain shape on the 9.7 fj.m silicate profile. Comparing this to 
the observations (Fig. [TTT i. the observed variations in molecu- 
lar sightlines might be explained if the dust grains in molecular 
clouds are more regular than in the diffuse ISM. 



Fig. 18. upper: The effect of changing the grain shape on the the 
T9 7 versus ^(J-K) relationship. The shape of the particles is rep- 
resented by a distribution of hollow spheres (DHS), where f„,ax 
represents the size of the cavities in the grains and is a measure 
for the irregularity of the particles (fma.x-^ means highly irreg- 
ular particles, f„,ax-Q. \ means al r nost sp herical particles). The 
Galactic Centre fit bv iMin et aP (l2007l) uses fma.v— 0.7. lower: 
The effect changing the maximum grain shape on the 9.7 //m sil- 
icate profile. The left panel shows the calculated 9.7 fm\ silicate 
absorption profiles in terms of the mass extinction coefficient 
and the right panel shows the same profiles normalised to 1 at 
10.5 jum as was done with the observations. The solid line in 
both panels represent the best fit for the Galactic Centre feature 
from lMinetal](l2007h . 



a mi x ture of amorphous 
T97I iGillett et all I l975t 



4.3. Chemical composition 

Interstellar silicate dust is mosth 
oUvines and p yroxen e s (e.g. iDav 
iKemper et all l2004t iMin et all 
Figure [19] shows the effect of changing the chemical composi- 
tion of the silicates on the T9 7 versus ^(J-K) relationship. This 
relationsh i p is pl otted for the chemical composition of the best fit 
iMin et all (l2007h found for the Galactic Centre, where we used 
graphite instead of amorphous carbon, at a mass fraction of 3 1 % 
(solid line; Draine, 2003). The dashed and dotted lines represent 
models with the same amount of graph ite, together with a type 
of amorphous pyroxene (i.e. MgSiOs; iDorschner et al.L Il995h 
only; or graphite and a type of amorphous olivine (i.e. Mg2Si04; 
iHenning & St ognienko. 1996) only, respectively. Going from a 
pyroxene to olivine composition rg 7/£'(J-K) decreases, but even 
assuming a pure olivinic composition is not sufficient to explain 
the observations. The best fit to the diffuse ISM feature towards 



the Galactic Centre contains Mg-rich end-members of silicates 
only (Min et al., 2007), and the contribution to the NIR opacity 
of these silicates is negligible. The presence of iron-containing 
silicates, or dirty silicates, would increase the contribution to the 
NIR opacity due to silicates, and thus decrease the amount of 
graphite required. 

Figure[T9](lower panel) shows the effect of changing the sili- 
cate composition on the shape of the 9.7 pm silicate feature. The 
change in shape of the 9.7 pm silicate feature that we observe 
could be explained if the silicate grains in molecular clouds con- 
tain relatively more amorphous pyroxene. This would increase 
the Tg 7 versus ^(J-K) relationship slightly, in contrast to the 
observations, but, since the observed variations in the 9.7 pm 
silicate profile are small, the effect on the T9.7 versus ^(J-K) 
relationship is probably also small. 
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Fig. 19. upper: The rg 7 versus ^(J -K) relationsh ip for the 
chemical composition of the best fit Min et al. ( 2007) found for 
the Galactic Centre together with a mass fraction of 31 % of 
graphite (solid line). The dashed and dotted lines represent a 
model with the same amount of graphite together with a pure 
magnesium rich amorphous pyroxene (i.e. MgSi03) and magne- 
sium rich amorphous olivine (i.e. Mg2Si04), respectively, lower: 
The 9.7 yum silicate absorption profiles for the above models. 
The left panel shows the calculated 9.7 yum silicate absorption 
profiles in terms of the mass extinction coefficient and the right 
panel shows the same profiles normalised to 1 at 10.5 jum as was 
done with the observations. The solid line in both pa nels repre- 
sent th e best fit for the Galactic Centre feature from iMin et all 
(I2007h . 



since above this grain size the 9.7 yum silicate profile changes in 
a way that is inconsistent with the observations. 

In the above models, we have not considered chemically in- 
homogeneous grains. For instance, metal inclusions can have an 
effect on the silicate absorption. Such changes should specif- 
ically affect the silicate optical depth, as the NIR extinction 
law remains un changed over a large range of column densities 
(iRoman-Zuniga et al.^.2007i) . 

5.2. E(J-K) as cause of the variations 

Since we did not find a model that can explain the breakdown of 
the T9 7 versus ^(J-K) relationship in terms of variations in T9.7, 
we consider it more likely that this is due to changes in ^(J-K). 
There are several ways to achieve such variations without af- 
fecting the shape and strength of the 9.7 fim silicate absorption 
band. 

The effect of grain growth and shape changes in graphite 
grains on the T97 versus ^(J-K) relationship has been investi- 
gated in Sect. ID Grain growth can be the answer ;/ silicates do 
not grow beyond about 1 fim and at the same time this growth 
should increase the ^(J-K) by up to a factor of 4. It is unclear 
what could be the cause for such a large increase in ^(J-K) in 
combination with modest silicate grain growth. Possibly, inho- 
mogeneous aggregates could cause such behaviour. 

In our model, the graphites and the silicates are considered 
to be separate grains, but in reality they might form composite 
grains. Moreover, our model does not include possible metallic 
iron grains, which could cause extinction in the near-infrared. 

The ^(J-K) in molecular clouds could also increase by 
the formation of a dust species that is less abundant in the 
diffuse ISM and has a high opacity at near-infrared wave- 
lengths, e.g. FeS, whi ch is thought to form in molecular clouds 
(IScappinietal.Ll2003h . 



6. Conclusions 

In order to study the dust properties in different environments in 
the ISM, spectra were selected from the Spitzer archive of highly 
reddened stars in both diffuse and molecular sightlines. We have 
analysed the 9.7 and 18 fim silicate absorption features in these 
spectra separately. Our main results are: 



5. Discussion 

5.1. Silicates (t^j) as cause of tiie variations 

The observed breakdown of the Tg 7 versus ^(J-K) relationship 
in molecular clouds could either be caused by a decrease of T9.7 
or an increase of ^(J-K). In order to explain the behaviour of 
the molecular cloud sightlines in terms of variations in T9 7, the 
silicate optical depth must decrease by at least a factor of 2 com- 
pared to the diffuse ISM silicates. Our model calculations show 
that such a large change in T9 7 is hard to achieve without sub- 
stantial changes in the shape of the 9.7 //m silicate absorption 
band. The observed modest band shape variations in molecu- 
lar cloud sightlines suggest that variations in rg 7 can not be the 
dominant cause for the difference between diffuse and molecular 
sightlines. Instead, variations in ^(J-K) may be responsible. 

The observed small changes in the silicate band shape are 
probably not caused by grain growth, unless grain growth causes 
more spherical grain shapes (see Sect.|4]i. Note that grain growth 
can still take place, but it is limited to a grain size of about 1 /im. 



- For diffuse lines-of-sight the strength of the 9.7 ixm silicate 
feature represented by the optical depth at about 9.7 yum (T9.7) 
shows a tight linear correlation with the near infrared colour 
excess, ^(J-K). Ho wever, this c oiTelation breaks down for 
molecular sightlines (IChiar et al.L [2007) . The measurements 
for the spectra analysed in this paper are consistent with this 
finding. 

- The shape of the 9.7 yum silicate feature is very similar for the 
three diffuse sightlines investigated in this study. Their shape 
is, within errors, also similar to th at observed in the lin e-of- 
sight towards the Galactic Centre (iKemper et al.L l2004l) . We 
conclude that we only observe small variations less than 0. 15 
in normalised t in the 9.7 yum silicate band shape in diffuse 
lines-of-sight, and that the Galactic Centre line-of-sight is 
representative for the diffuse ISM in general. The shape of 
the 18 yum silicate absorption feature is also, within errors, 
found to be invariable for the diffuse sightlines. The situa- 
tion is different for lines-of-sight passing through molecular 
clouds: we observe small variations on the short wavelength 
side of the 9.7 jim silicate band. 
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- Since both the shape of the 9.7 /zm silicate feature and the 
relationship between T9.7 and ^(J-K) change for molecular 
sightlines as compared to diffuse sightUnes, we conclude that 
the dust properties in molecular clouds are different from the 
diffuse ISM. 

- Based on comparison with theoretically calculated profiles, 
we have ruled out grain growth beyond about 1 fim as be- 
ing responsible for the changes in the 9.7 //m profile. The 
observations could be explained if the dust grains are more 
spherical in molecular clouds than in the diffuse ISM, but the 
underlying process is not understood. 

- We propose that the flattening of the T9.7 versus E{J-K) re- 
lationship is caused by an increase in E(}-K) in molecular 
clouds. The mechanism that causes this, however, is not yet 
clear. Further research is necessary to investigate whether it 
can be caused by grain growth. 
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